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Abstract 

We introduce a non-overlapping variant of the Schwarz waveform relaxation algorithm for semi- 
linear wave propagation in one dimension. Using the theory of absorbing boundary conditions, we 
derive a new nonlinear algorithm. We show that the algorithm is well-posed and we prove its conver- 
gence by energy estimates and a Galerkin method. We then introduce an explicit scheme. We prove 
the convergence of the discrete algorithm with suitable assumptions on the nonlinearity. We finally 
illustrate our analysis with numerical experiments. 

1 Introduction 

Schwarz waveform relaxation is a new class of algorithms for domain decomposition in the frame of 
time dependant partial differential equations. They are well-adapted to evolution problems, designed 
to solve the equations separately on each spatial subdomain on the whole time interval , exchanging 
informations on the space-time boundary of the subdomains, overlapping or not [BJ. In particular for 
wave equations, it is of great importance, due to numerical dispersion, to be able to handle local time and 
space meshes, and this is allowed by the present method. We presented the method for the linear wave 
equation in 7\ and [5] . When using overlapping subdomains and "classical" Schwarz waveform relaxation 
-by a Dirichlet exchange of informations on the boundary- the so defined algorithm converges in a finite 
number of iterations, inversely proportional to the size of the overlap, which can be penalizing. We 
introduced optimized transmission conditions, relying on the theory of absorbing boundary conditions, 
which improve drastically the convergence of the algorithm. 

Very little has been done so far about nonlinear Schwarz algorithms. An analysis of the classical 
Schwarz waveform relaxation algorithm was performed in [4] for a conservation law. The goal of this paper 
is to define new Schwarz waveform relaxation algorithms for the semilinear wave equation. We introduce 
two nonoverlapping algorithms. The first one referred to as linear uses the absorbing boundary condition 
of the linear problem whereas the second one referred to as nonlinear uses the nonlinear absorbing 
boundary conditions designed by J. Szeftel in [TT] . 

In Section 2, we introduce the definitions of the algorithms. 

In Section 3 , we prove the algorithms to be well-posed. For the precise analysis, we use a fixed point 
algorithm with regularity estimates on a linear problem. 

In Section 4, we prove the convergence of the algorithms. The proof is an extension of a clever trick in 
.8], already used for linear algorithms, either hyperbolic or parabolic (see [7])- However the nonlinearity 
requires a very fine analysis. 

In Section 5, we design discrete Schwarz waveform relaxation algorithms. In each subdomain, the 
interior scheme is the usual leapfrog scheme for the linear part, with a downwinding in time for the 
nonlinear part. The exchange of informations on the boundary is naturally taken into account by a finite 
volume strategy. In Section 6 we study the convergence of the algorithms, by discrete energy estimates. 
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As it is always the case for nonlinear problems, the well-poscdncss and convergence results hold only 
locally in time. Therefore numerical experiments are very important to bypass the limitations of the 
theory. We present the results in Section 7, showing in particular that our nonlinear algorithm gives 
optimal results within a large class of algorithms. 

Remark Due to the complexity of the mathematical theory, we restrain ourselves to the one dimensional 
case. The multidimensional study contains additional difficulties due to the geometry and should be the 
heart of a forthcoming paper. 

2 Problem Description 

We consider the second order semilinear wave equation in one dimension, 

{d 2 t -d 2 x )U = f(u,dtU,d x U) (1) 
on the domain R x (0, T) with initial conditions U (•, 0) = p, dt U (•, 0) = q. 



2.1 Absorbing Boundary Conditions for the Semilinear Wave Equation 

The question of absorbing boundary conditions arises when one wants to make computations on an 
unbounded domain: a bounded computational domain is introduced, on the boundary of which boundary 
conditions must be prescribed. These boundary conditions must be absorbing to the waves leaving the 
domain. A whole strategy has been designed by Engquist and Majda for linear problems with variable 
coefficients, using pseudo-differential operators [3J. Recently it has been extended to nonlinear operators 
by J. Szeftel, in particular for the semilinear wave equation using the paradifferential calculus of [lj 
and [9] . We introduce a family of operators 

B ± (g ± )u = d t u±d x u + g ± (u), (2) 

for C°° functions g^ such that <? ± (0) = 0. The linear absorbing boundary operators are given by g^ = 0. 
In the case where f(u,Ut,u x ) — fi(u) + f2(u)u t + fz{u)u x with fj in C°°(R), 1 < j < 3, and /i(0) = 0, 
the following nonlinear boundary operators are given in [llj : 

i r u i r u 

9 + {u) := -- J {h / 3 )(fK. 9~{u) ■= -- j {h + h)(m- (3) 
We replace the problem on the domain R by a boundary value problem in f^o = (a, b): 

(df - dl)u = f(u, d t u, d x u) in n x (o, t), 

B-{g-)u(a,-) = 0, B+{g+)u(b, ■) = 0, U 

with initial values p and q. Such boundary conditions give well-posed initial boundary value problems, 
and are absorbing provided the intial data be compactly supported in ^o, see Following the strategy 
in [7], we use such absorbing operators for domain decomposition. 

2.2 A General Non-Overlapping Schwarz Waveform Relaxation Algorithm 

We decompose the domain (a, b) into / non overlapping subdomains Qi = (a^, a i+1 ), aj < ai for j < i and 
ai = a, = b, and we introduce a general non overlapping Schwarz waveform relaxation algorithm. 
An initial guess {/i^'°}i<j<i+i is given. For k > 1, one step of the algorithm is 

" (d 2 t - dl)u\ = f(utd t u1,d x u1) in fit x (0,T), 

uf(-,0)=p, a t uj(-,0) = ginn <) 
_B-(g-)uHa i ,-)=h7< k -\ B+(g + )u>?(a i+ i,-) = hi' k - 1 m(0,T), (5) 

hr- k =B-{g-)u k _Mr), h+> k = B+(g+)u k +1 (a i+1 ,-) in (0,T), 
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where B^ are given in ([2]). For ease of notations, we defined here hf' k = and hjl_i = 0, so that the 
index z in ([5j) ranges from i = 1, 2, . . . , /. In the sequel, we call linear transmission condition the choice 
g ± = and nonlinear transmission condition the choice ([3]). For the classical linear homogeneous wave 
equation, it has been proved in [7] that the algorithm converges optimally if T is small enough (which 
means in two iterations, independently of the number of subdomains), and the transmission operators 
Z?* are given by S ± = dt ± d x . This behavior is due to the finite speed of propagation, together with 
the fact that these operators are the exact Dirichlet Neumann operators in this case. In the nonlinear 
case, the propagation still takes place with the finite speed, but we can use only approximate Dirichlet 
Neumann operators. Therefore the classical Schwarz algorithm with overlap is still convergent, and for 
our nonoverlapping nonlinear algorithms, we will use energy estimates. 



3 Well-posedness For The Subproblems 

The study of the nonlinear problem relies on an iterative linear scheme. Therefore a first step for the 
definition of the algorithm is the study of the nonhomogeneous initial boundary value problem for a 
general domain SI — (a_, a + ), 

(Of - 8 2 x )u = f(u, d t u, d x u) in fi x (0, T), 

B-(g-)u(a-,-) = h-, B+(g+)u(a+,-) = h+, 

with initial values p and q. We will use for j ' < 2 the spaces 

Vj(Q, T) = {ue L°°(0, T; L 2 (Sl)), d a u € L°°(0, T; L 2 (Sl)), \a\ < j}. (7) 

In formula 10, a is a 2-index in N 2 , the first coordinate in a stands for the time, and the second one 
stands for the space, so for instance d a u = dt x u for a — (1,1). Vj(Cl,T) is equipped with the norm 

IM[yj(n,T) = m.ax\ a \<j ||<9 a u|| L oo (0:T;L 2 (f2 )). 

Theorem 3.1 Let p in H 2 (fl) and q in H (SI). There exists a time T* such that for any T < T* , for 
h in i2 rl (0,T) with the compatibility conditions 

h ± (0) = q(a±)±p'(a±) + g ± (p(a ± )) 1 (8) 

(0|) has a unique solution u in V2(Sl,T), with dtu(a±,-) and d x u(a±,-) in _ff 1 (0,T). Furthermore there 
exists a positive real number C* such that 

NIWt)+E E \\d a u(a ± ,-)\\ 2 W) <C*(|H|^ (n) + ||g||^ 1(n) 11^11^(0,^), (9) 

± |q| = 1 ± 

where T* and C* depend on the data p 1 q, f, g^ . 

This result has been first proved in [11] with homogeneous boundary conditions (i.e. h — 0). The 
additional difficulty comes from the boundary conditions, and we give here the main steps of the proof. 
It relies on the construction of a sequence of linear problems of the form: 

d?u-dlu + u = F in Six (0,T), , . 

(d t u-d x u)(a.,-)=H-, (d t u + d x u)(a+,-)=H+. [U) 



Proposition 3.2 Let p in H 2 (S1) and q in H 1 (Sl). For any positive time T, let F in i? 1 ((0, T) x SI), 
and in iJ^OjT) with the compatibility conditions 

H ± (0)=q(a ± )±p'(a ± ). (11) 

Then, &10\) with initial data p and q has a unique solution u in V2(Sl, T), with dtit(a±, •) and d x u(a±, •) 
in -ff 1 (0, T). Moreover we have the following bounds on the solution 

ll"lly 2 (o,T) +J2(\\ d Ma±,-)\\ 2 H^o,T) + \\ d xu( a ±r)\\ 2 Hi(o,T)) 
± 

<Cie T (ll^ll^(o,T^(n)) + ll^(-.0)||l 2(n) +^11^11 , (12) 

± 

where C\ is a universal constant. 
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Proof We start with the a priori estimates. We multiply (|10|) by dtu and integrate by parts in 
1 d 



2 dt 



{\\dtu{-Mh(n) + \H-Mm ( n))+{d t u{a-,t)) 2 + {dtu{a + ,t)f 



[F(; t), d t u(-,t)) L 2 m + H-(t)d t u(a-,t) + H+(t)d t u(a+,t). 



Using the Cauchy-Schwarz inequality on the right hand side, together with the inequality a/3 < ^ a 2 + ^ /3 2 
for all a, [3 £ R, and finally integrating in time, we obtain 



\\dtu(;t)\\i H n) + \H;t)\\H H n } + / [(d t u(a-, s)) 2 + (d t u(a+, s)) 2 } ds 

Jo 

< f\\\d t u(; S )\\ 2 L2(n) d S + f \\F(; S )\\l 2{n) ds 

Jo Jo 

+ \\Q\\l H n) + \\p\\m ( n) + / W(*)) 2 + (H + ( S )) 2 } ds. 

Jo 

By Gronwall Lemma, we deduce that 



\\dtH-,t)\\h m + \\H-,tW HHn) + / [{d t u(a^,s)y + (d t u(a + ,s)Y}ds 

Jo 

+ \\q\\h(n) + \\p\\m(n) + W H± Wh(o,T)), 

± 

which gives 

max ||^ a u||i ! x.(o,T;L2(o )) + \\dtu(a^,-)\\ 2 L2{Q T) + ||<9 t u(a+, •)|| 2 2(0 T) 

< e T (ll^ll^ ((0 ,T)xo) + hWhm + \\p\\m { n) + £ H^lliW" (13) 

± 

Differentiating in time in (|10j) . we now apply (|13p to dtu, and obtain: 

max \\d a d t u\\ 2 L ^ { Q tT . L 2 [n)) + \\d 2 u(a-, -)lll=(o,T) + \\ d tu{a+, OIll^o.T) 

< e T (\\d t F\\h {i o,T)xn) + \\dM-M\h { n) + h\\m(n) + E W d * H± I^W" (14) 

± 

We must estimate dfu(-,0) in the righthand side of (fT4|) . We multiply (110|) by dfu, integrate in space 
and evaluate at time 0: 

\\dM-M\h m + (d xP ,d x d 2 u(;0)) + (p,^fi(-,0)) 

+ q(a-)d 2 u(a^, 0) + q(a + )d 2 u(a+, 0) 

= H-(0)d 2 u(a.,0) + H+(0)d 2 u(a + ,0) + (F(;0),d 2 u(-,0)). 

We integrate by parts in the second term, and rewrite the equality as 

\\dM; 0)||| a(n) - (d 2 u(; 0), d 2 u(-, 0)) + (p, d 2 u(., 0)) 
= (F-(0) - ?(o_) + 9 a p(a_))0 t 2 u(a_,O) + (#+(0) - g (o+) - d xP (a + ))d 2 u(a + , 0) + 0), c\ 2 u(-, 0)). 

The boundary terms on the right-hand side vanish by the compatibility conditions, and we get 

Hu{-M\l* m = (d 2 x p-p + F(;0),d 2 u(-M- 



4 



Using the Cauchy-Schwarz Lemma, we obtain 

\\d*u(;0)h*m ^ \\%P-P + n,0)\\»piy 
We replace the term \\d^u(-, 0) ||_L 2 (n) m (HHJ) , and we deduce the second a priori estimate: 

max \\d a dtu\\ 2 LX{Q T . L2m + ||<9 t 2 w(a_, -)lli2 (0iT) + ||<9 t 2 u(a + , Olli^o.r) 

<3e T (||c> t F||| 

2 ((o,t)xo) + II -^(^ 0) lli,2( ) + ||p||ff2(si) + 11911^1(0) + E H 5 *- ff± lli 2 (o,T))- ( 15 ) 

± 

We still need to estimate the mixed derivatives d xx u in the interior and d x tu on the boundaries. We use 
the equation, which gives in the interior 

||0ra^||L°°(O,T;L 2 (n)) < \\dttu\\ (0 ,T :L 2 (Q.)) + INU°°(0, T;L 2 (fl)) + \\F\\ L°°(0,T;L 2 (fi)) • 

We now introduce the inequality 

H- F llz,<-(0,T;L 2 (O)) < 2 (ll- F (->0)lli 2 (n) + T ll 5 t^llL 2 (0,T;L 2 (n))): 

and by (p~5|) and (fl3|) . we get, using that e T > 1 and e T > T, 

||9 2:a; M|| i oo( 0iT;L 2 (n) ) 

2 H 1 (0,T))' 

± 

As for the boundary term, we get for instance on the left boundary 

\\d tx u{a„,-)\\ L 2 {0 T} < \\d tt u(a-, -)IIl 2 (o,t) + \\d t H^\\ L 2 {0 T} . 
Squaring the inequality, and adding the term coming from the right boundary leads to 

Ell 9 *-"( a ±'-)lli 2 (o,T) <2^(||5 tt? i(a ± ,-)|li2 (0iT) + ll^ll^^)), 
± ± 

which provides the last estimate announced in the proposition. The well-posedness is then derived in a 
standard way by the Galerkin method. ■ 
The solution u of the nonlinear subdomain problem is now defined through an iterative scheme. The 
initial guess is uq — p. At step fc, Uk being known, we define 

F(w) = w + f(w,d t w,d x w), g ± (w)=g ± {w(a ± ,-)), H± (w) = h ± - Q± (w) . (16) 

Uk+i is the solution of the linear initial boundary value problem (jlOp with data = J-{uk), g^ = G^iuk), 
= 7i(ufc), and initial data p and q. The proof of convergence for the sequence is written in details 
in [11) . The uniqueness follows from the result: 

Lemma 3.3 There exists a real positive increasing function 9 such that, for any time T , for any v in 
V 2 (fL,T), d a T{v) is inVi(n,T), Q ± {v) andH{v) are in L°°(0,T). Moreover, for v lt v 2 inV 2 (tt,T), we 
have 

H^Oi) - U ± {v 2 )\\loo{o,t) < 9{\\vi\\v 2 (n,T) + ||«2||v 2 (n,T))IK - v 2 \\ Vl (n,T), 

(17) 

W^ivi) - F(v 2 )\\ Vl ( Q ,T) < 9{\\vi\\v 2 (n,T) + IKIk 2 (n,T))IK - v 2 \\v 2 (n,T), 
As a consequence, we have the well-posedness of problem 

Corollary 3.4 Let p in Hq(£1) and q in ffg(O). There exists a time Tq such that for any T < Tq , 
has a unique solution u in V 2 {Q,,T), with d a u(a,-) and d a u(b,-) in H l (0,T) for \a\ — 1. Furthermore 
there exists a positive real number C* depending only on the size of f2 such that 

NW)+ E \\^M\\h(o,T) E ll^.-)ll^(o,T)<C*(IIPllH 2 (a) + lkll^ ( n))- 

\a\ = l |a| = l 
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4 Convergence of The Algorithm 

We now study the convergence of the Schwarz waveform Relaxation Algorithm ([5]). In order to define 
the algorithm, we need a regularity result: 

Proposition 4.1 For any e, < e < 1, V 2 (Q,T) C C°(0, T; i/ 2_£ (Sl)) n C^O, T; iJ 1_e (f2)). 

Proof By using extension operators in time and space, it suffices to prove the result in I x 1. We make 
use of the Littlewood-Paley theory (see for example [2]). In particular, there exists tp and \ two tempered 
distributions on R, with ip supported in (—8/3, —3/4) U (3/4,8/3), x supported in (—4/3,4/3), and 

x(£) + 5>( 2 ^) = 1 > v £ elEL 

We define the dyadic projectors A 9 by their action on a function u, 

A_iw = x(D)u, A q u = tp(2- q D)u for q > 0, (18) 
where D — —id. These operators give an equivalent norm in .fP(R), 



E 2 2?S ||A ( 
?>-i 



u\\ 2 



They can also be used to define the Zygmund spaces 

CI = {ue S', I u| r = sup 2 <?r ||A, 3 u|| 2 < +oo}. 
?>-i 

CI coincides with the usual Holder space when r is not an integer. For any positive r, we know that 
W r, °°, the space of functions in L°° with derivatives of order up to r in L°°, is included in C£. Therefore 
we have 

y 2 (M,R) C C°(R,ff 2 (R)) nCjftff^R)) nC 2 (R,L 2 (R)). 
We need an interpolation lemma. 

Lemma 4.2 For any positive a, (3, a, b, for any 9, < 9 < 1, 

C(R,# Q (R)) nCf(R,ff b (R)) c C» a+(1 ~ 9)/3 (R,tf 9Q+(1 - e)b (R)). 

Applying the lemma with successively (a, (3, a, b) = (0,1, 2, 1) and (a, /3, a, 6) = (1,2,1,0), we find for 
any 9,9' in (0,1), 

y 2 (R,M) c C, 1_e (R,iJ 1+e (R)) nC 2 ~ e '(R,i? e '(R)). 
Since for any e > we have C| C C° and Cl +e C C 1 , this concludes the proof of Proposition 14. II 

Proof of Lemma 14.21 It relies on the convexity of the exponential function. 



|,.||2 _ n 2 2j(fla+(l-e)0) 2 2/c(fla+(l-fl)b|| A t ax 



J^- 1 fc>-l 

where A*- (resp. A|) is the Littlewood Paley operator acting in the time (resp. space) variable. 



>z„.l|2 



fe>-l fc>-l 

< 



E 2 2fc ( ea+ ( 1 - e ) f, ||A*A>|| 2 = E ( 22fca H A l A >H 2 ) (2 2fcb ||A*A 

" :-l V 

E 2 2fca ||A*A^|| 2 )7E 2 2fcb ||A5A^|| 



fc>-i fe>-i 
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Therefore we have 



IMI^ + a-^ (K ,H 9 « +(1 -*>» (R)) < (<™P >/ ja E 2 2fc 1|A*A^f) e (sup 2^ £ 2 2 ^|| AjA^I 2 ) " 
which writes 

||u||^ Q + (1 _ W(R ^ o+(1 _ e)()(R)) < (|Mlc ? (R„ff»(R))) (iNj 



Theorem 4.3 Le£p m Hq(£1) and q in Hq(H,). There exists a time T\ < Tq such that for any T < T\, for 
any initial guess hf in H 1 ^, T) with the compatibility conditions hf(0) = q(aj+i) +p' (aj+i) + g + (p(aj + i)) 
and ^7(0) = li a i) ~P'( a i) + 9~{p( a i))i th e algorithm is defined and converges in UiV2(£l;,T) to i/ie 
solution u of 0). 

Proof We first prove that the algorithm is well-defined : with the assumptions on hf in the theorem, 
we know by Theorem 13.11 that (0 defines in each fli a u\ in V^(fij,T), with dtu\{di,-), dtu\ (flj+i, •), 
d x u}(ai, •) and ftu^tii+i, •) in ff^O, T) for T < X*. Furthermore, by Lemma I3T31 B - (5 _ )r4_i(ai, ■) and 
B + {9 + ) u i+i( a i+ii ') are m Lf 1 (0,T). As for the compatibility conditions, we have 

S _ (p _ )iij-_ 1 (a f) 0) = Inn (ft u^a^i) - x «i_i(a<,*) 
and by Proposition [4TT] we can pass to the limit and get 

# _ (ff _ )«i-l( a ii°) = g( ffli ) -p'(Oi) + S~(K a i))- 
This, together with the same regularity result on aj+i, permits the recursion. 

We define for T < min(Tg , mini(T* fc )), for k > 1, the quantities (with u% — u/nj for 1 < i,j < I, 
j = i or j = i - 1, 

e* = u * - m, 
f t k = T{ut) - H^i) 
h L 3 = 9~(Uj(ai)) - h^f = g+(u k {a i+1 )) - g+ ( Uj (a l+1 )). 

The operators T is defined in (TIT))) . The error e\ in fi^ at iteration k is a solution of 

(ft 2 -ft 2 )ef + e k i= ff mCliX(0,T) (19) 
(ft - ft)e? + ftfr" = (ft - ft^ 1 + <7_V on {a t } x (0, T) (20) 
(ft + ft) + = (ft + ft^ 1 + <7+i + on x (0, T) (21) 

with vanishing initial values and e$ = 0, ef +2 = 0, /i^' = 0, h^'^ +1 = . In order to get a new energy 
estimate in fii, we multiply (|19p by dte\ and integrate by parts: 

f t EnM) ~ [d t e%eUa i+ i, ■) - d t e%e^( ai , •)] = (/f,fte 8 fc ) (22) 

with Eq(u) — 5(||ftw||| 2 (f2) + llft^lli^n) + II u IIl2(o))- We rewrite the boundary terms using the boundary 
operators: 



ftefftef (a i+1 , •) = -((ft + d x )e*(a i+1 , •) + h^f- -((ft - ft)ef (a i+1 , •) + Z^) 2 + J$, +1) 

1 i 4 (23) 

-fte*ft4(a i( •) = ^((ft - ft)ef (a 4 , •) + /^"^((ft + d x )e\ (a*, •) + fe*L+ (i ) a + R* ti _ v 
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The remainders R k i+1 and R k i _ 1 will be evaluated later. We insert (|23|) into (|22|) , and obtain 

Using the transmission conditions (f2T)|) , (f2T|) . we get 

J t E ^) + ~ d x )e k M+x,-) + ^-Im) 2 + + d x )e k M, + 

(24) 

We sum ([24]) on the indexes i, 1 < i < I and integrate in time. We translate the domain indexes in the 
right-hand side. Defining 

EonAe k ) = \m - d x )eHa i+ i, •) + h k ^f + ((d t + d^fa ■) + /£+ 4 ) 2 ], 
we get, since the initial data vanish, 



i=i - 70 i=i - 70 i=i ^ 

+ / E(^+i + ^-i)( s ) ds + / Eat^n^. 

Differentiating the equation and the transmission conditions in time yields the bound on dte k : 

i=i - 70 i=i - 70 i=i 

+ / E(^t+i+^i-i)( s ) rfs + f i^idtf^^due^ds. 
Jo l=1 Jo t=1 



(25) 



(26) 



We now estimate the remainders. We start with Rf i+ i (ignoring the superscript k): 

= 2i- h tr( 2 ( d t e t + d x ei)(ai+u + + hr +l l (2(d t ei - d x ei)(a i+1 , •) + V +M )], 
and we get a bound on the integral of Ri^+i. 

We can treat Ri^-i,Ri^-\, and Ri y i-\ the same way and obtain 

iZ,,i-i(s)ds < |(|| Vifi +||^tliH 2 I )+ 1 1 {\\d t e l {a tr )\\ 2 x + Re l (a l ,-)|| 2 i \ 

4 M "H5(0,t) ' ffl(0,t) 2 H 2(0, t) "'H 2(0, t) 

Ai, m («)d«< KWdthfi || 2 , +||9 t V +1 J 2 , ) + i(||a 2 ei (a l+1 ,-)|| 2 . +11^(^+1, -)H 2 i ) 
Xi-i(«)^< KH^Vill 2 i, +l|3At li || 2 . ) + i(||a 2 ei ( ai ,-)|| 2 > + ||^ t e,(a,, -)ll 2 i, 



At point a,i for instance, by the Trace Theorem, there is a constant C3 independent of T, such that for 
any a with \a\ = 1, we have 

ll^cOI^-i^ < ||fl 0, e i (o < ,-)||L«(o, t )<C3||fl 0, e < || i ri(n 4X (o ) t)), 

\\ dad M^-)\\ H - im < \\Vei(fi<>-)\\ H i m <C 3 ||9 Q ei |l^(n,x(o, t )), 

which gives our first bounds on the remainders: 

f R t , l+1 (s)ds < l(\\ht,if Him + W h i+i,i\\ 2 Him ) + ^INI^(n ( x(o,t)), 

Ri,i + i(s)ds < l(\\d t hf. \\ 2 h i ^ + \\dtK +1 .\\ 2 H i m ) + ^-WetWl^^y 
We now insert the previous estimates in ([25]) and (T3H|)- By Cauchy-Schwarz inequality we get 

1 f t 1 

$>n,(e 4 fc )(i)+ / J2 E ™Aei)(s)ds< 

J ° 1=1 8=1 1=1 ^° 

i — 1 i — 1 

Jo l=1 

E^^^^+^^Eii^ii^^x^ + ^E f Hfi\\h m (s)ds 

i=l i=l t=l ^° 



< 

i=l i=l 

+ 1 E(i^^iik ( o, t) + n^AW + 3 E(ii^n^ (M + \\^r Bi ) 

i—1 1 

Adding (gTJ) and (J2SJ), we can write 

£(£ ni (e*) + E ni {d t e\)){t) + / ^(£; a0l (e 4 fc ) + WftefrJOOda 

i=l ■'0 i=l 



< /' E(^ ( e " _1 ) + ^ (fte*- 1 ))^)^ + (Cf + 1) I 
" , ° i=i i=i 



e i\\ 2 H 2 (niX{0,t)) 



3 1 



- 2 E(Kf i^ (0 , t) + ifeu^ (0 , t) + Kn^ (0)t) + n^t«ii^ (0 , t) ) 

2—1 

3 J 

-4 Z^"""'* \\ H i m t ll^+i,; ii^j -r nuii^ n H i (o t) -r ll"*»i-i,ill ff l (0it) J 

2—1 



+ IIWI L 2 (f!x(0,T))i 



i=l 



(28) 



We now estimate the quantities involving the h*'^. Refining the results in Lemma 13 . 31 we have a real 
positive increasing function 62, such that 

ll^fH^^ + H^Ukco.t)^^^ E (ll^ll^(^x(0,T)) + l|9X fc Hi 2 (nx(0,T)))) E H 9Qe "H^(Ox(0,T)) 

|c|<2 | a |<2 
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which gives 
I 



Y^iEuM) + EnM^W) + I ^2(E ani (e^+E ani (d t e^)(s)ds 
i=i Jo (=1 

< fjZ^^ 1 ) +E d n i (d t et 1 ))(s)ds 
J o i= i 

7 

+ E 6 ' 3 ( X! (ll 9O<U *lli 2 (fix(0,T)) + H 5Q "! i: |li 2 (Ox(0,T))))ll e 4 fc |lff2(o I x(0,t)) 
i=l |q|<2 
<*#(0,2) 

L 2 (nx(0,T)) 

with 3 = 36*2 + C| + 1. We now evaluate the terms in the right-hand side. We first note that 
IWHh^m) < /W(e?) +Eo i (8 t e$) + ||S»xe?|[£ a(n0 )ds, 



and evaluate \\d X xef\\ 2 L 2 (n . y(Qt )) by equation (JTjJ 



i llL 2 (niX(0,t)) 

ll^aKc e » llla(n iX (o,t)) - 3(||t? ti ef 11^2^. x (o,t)) + \\ e i\\l 2 (n z x(o,t)) + l/* fc lli 2 (o, x (o,t)))' 
from which we deduce 



rt 1 



(29) 



L 2 (OiX(0,t))- 

There remains only in (|29p 

J2(EnM)+E Qi (d t e^)(t)+ / J2(E aQi (e?) + E aQi (fte? ))(s)ds 
i=i J ° i=i 

,t £ 

l 2 (Qx(0, 

1=1 |a|<2 *' 

o^(0,2) 

L 2 (Ox(0,T)) ' II "i lli 2 (nx(0,T)) )) + 2)dl/iH L 2 (f2x(CLT)) + ll^t/j lli 2 (nx(0,T)))- 

1=1 M<2 

a^(0,2) 

(30) 

Again, as in Lemma l3.3( there exists a positive increasing function #4 such that 

011^(0) + Hfi(t, Olll^n) < ol( E (H^lli 2 (o) + liaVlliW E H^illW 

c|<2 | Q |<2 

«#(0,2) "#(0,2) 

the latter sum means that no term d xx are present, therefore we can bound the sum by twice the energy. 
Furthermore we know that the energy of u is bounded on the interval (0, T). Thus there exists a new 
positive increasing function #5, depending on u, such that 
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We insert (f31"j) into (|30p . and get (with a new function 9q ) 



^(^(e?)+^ t (fte?))(t)+ / £(ifcn«(e?)+2WfceJ))(s)ds 
i=i J ° i=i 

•'o i=l 

+ 2 / (0 6 (i?o I (e 4 fe )+^(a t e J fc )))(^(e J fc )+£;n I (a t e 4 fc ))( S )d S . 

For clarity we define 

EL = (e?) + ^ (fte*)), ^ = (e 4 fc ) + ^ (fte?)). 



i=l »=1 

and we can rewrite (1321) as 



(32) 



i& t (t)+ / ^'(s)d s < / £; 6 fe - 1 ( s )d S + / e 6 (El t ( s ))El t (s)ds (33) 

Jo Jo 
Summing in k, we define Ef nt = X)*?=i E\ ntl and we have 

E*t{t)+ ( E*(s)ds< [ E°(s)ds+ I 6 (EL(s))Ef nt (s)ds. 



Let now C > 0. If t tends to 0, J Q El(s)ds + tC6e(C) tends to 0. Therefore there exists a such that 
J Tl ££(s)ds + TiC6» 6 (C) = C, and so we have for t<T x , 

We conclude that u\ exists on the time interval (0,Ti), and that $Zf = i (e*) + E^dte^)) tends to 
when fc tends to infinity: the sequence uf converges to u on (0, T\) in each subdomain in the norm of 
energy. ■ 



5 A Finite Volume Discretization 

We use here a finite volumes scheme, which has been described in [7] for the linear one-dimensional wave 
equation, and extended to the non linear boundary value problems in the frame of absorbing boundary 
conditions in |10j . We restrict ourselves to uniform meshes in time and space. 

5.1 Discretization of the Subdomain Problem (EI) 

The domain CI x (0, T) is meshed by a rectangular grid, with uniform mesh sizes Ax and At. There are 
J + 1 points in space with Ax = (a+ — a_)/J, and N + 1 points in time, with At = T/N. We denote 
the numerical approximation to it(a_ + jAx,nAt) by U(J,n). We introduce the notations: 

Dt UU, n) = ^•,n + l) -l70-,n) [ ^ _ tffc „) - tffc „ - 1) 



At ' t w. y Af 

i* UU, n) = ^• + l,^-^,n) [ D _ ^ . n) = t /( J -,n)-^-l,n) | 

D o ^ B) = ^ + l,n) 2 -^0--l,n) > ^0 u{ . n) = UU,n + l)-UU,n-l) 

_ 3U(j,n)-4U(j,n-l) + U(j,n-2) f D~ U(j, n) for n > 2, 

^ U[j ' n) ~ 2At ' * U[J ' n) ~\ D^U(J, n) forn = l. 



(34) 
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The last finite derivative in (|34[) is a second order approximation of dtu, to be used in the nonlinear 
term, in order to design an explicit scheme. 
The scheme in the interior writes 

(D+ D- -D+D-) U(j, n) - f(U(j, n), £> t "* U(j, n),D° x U(j, n)) =0 in [1, J — 1| x [1, N - 1]. (35) 
We define the discrete initial value as 

P{j) =p{a- +j Ax), Q(j) = q(a_ +j Ax), 
and we obtain the initial scheme 

(D+ - D-)U(j,0) = Q(j) + ^-f(P(j),Q(j),D° x Q(j)), in|l,J-l|. (36) 

For the boundary conditions, we dehne the discrete boundary operators as: 

B~ (/, g- ) (7(0, n) = (D° -D+ + ^-D+ D~ )E7(0, n) - ^/(f/(0, n), D~* (7(0, n), D+ U(0, n)) + g~ (U(0, n)), 
B-(f,g-)U(0,0) = (D+ -D+ + ^L D +)U(0,0) - ^Q(O) - ^/(P(0), Q t (0), D+ P(0)) + g-(P(0)), 

B+U, 9 + ) U(J, n) = (D° +D~ + D~ )U(J, n) - ^f(U(J, n), D~* U t {J, n), D~ U(J, n)) + g+(U(J, n)). 

B+(f, g+) U(J, 0) = (D+ +D~ + )U(J, 0) - ^Q(J) ~ ^/(P(J), Qi(J), D~ P( J)) + g+(P(J)). 

We define the boundary data as 



(37) 



(38) 



^ ^t„+At/2 2 /•At/2 

H±{jl) = M h ± (r)dT, l< n <N,H ± (0) = — h ± (r)dT (39) 

The discretization of Problem © is now given by ([33]) , ([5^1 , with boundary conditions 

B-(f,g-)U(0,-)=H- , B+(f,g+)U(J,-)=H+ in [0,n</Y], (40) 
where the discrete boundary operators are given in (|37|) , (f38|) . 

Our numerical computations indicate that this scheme is second order both in space and time. 

5.2 The Discrete Schwarz Waveform Relaxation Algorithm 

The equation is now discretized on each subdomain fij x (0, T), i = 1, . . . , I separately, using an uniform 
mesh with sizes Ax and At. There are Ji + l points in space and N+l grid points in time in subdomain il^, 
with Ace = (ai + i — ai)/Ji and At — T/N. We denote the numerical approximation to itf (a^ + jAx, nAt) 
on £!i at iteration step k by U^(j, n). 

The problem in domain fli is now defined through boundary data iJ J ± ' fcl , coming from the neigh- 
boring subdomains f2i±i .Therefore, we define the extraction operator from domain fi, to his neighbours 



as: 



B+(f,g+)U t (0,n) = (D t ° + D+ - D~ )[/,(0, n) + ^/((7,(0, n), D~* (7,(0, n), D+ (7,(0, n)) + 9 +((7,(0, n)), 

B+(f,g+)U i (0,0) = (D+ +D+ - )l/,(0,0) + ^Qi(0) + ^/(P,(0), Q,(0), D+ P,(0)) + 9 +(P,(0)). 

(41) 

B-{f,g-)Ui{Ji,n) = (D° t -D- - ^D+D-)C/,(J„n) + ^/(C/,(J„n),Z3 t -*C/,(J„n),DJC7,(J„n)) + g-(C/,(J,,n)), 
B-(f,g~) Ui(Ji,0) = (D+ -D- - ^D+)Ui(Ji,0) + ^Qi(Ji) + ^f(Pi(Ji),Qi(Ji),D- P,( J,)) + g~ (P,( J,)). 



(42) 
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The discrete Schwarz waveform relaxation algorithm on subdomains fii, i = 1, . . . , / is defined as follows. 
An initial guess {H®' }i<;</ is given. For k > 1, we solve 

(D+ Dt -D+D-)Uf = f{Ut DT* Ut D° x U?) in [1, J t - 1] x [LA], 

U?(.,0) = P i; (A + - ^£ D- )U?(; 0) = Q 4 + ^f(Pi,Qi, D° x Qi) in [1, J, - 1], 

2 2 (43) 

B-(f,g-)Uf(0,-) = Hr^- 1 , B+ {f , g+)U>* {J h ■) = H^ 1 in [0,n< AT] 

Hr> k =S-(f,g-)U?_ 1 (Ji-i,-) , H+> k = B+(f,g+)U? +1 (0,-) in [0,n<A]. 

As in the continuous algorithm, we set Uq = and Uj +2 = 0. 

We denote by U the discrete approximation of problem (01, obtained by solving (|35l36l40p on Q = (a,b) 
with J = J% intervals of length Ax, and = 0. Each subproblem is an explicit scheme, thus has 

a unique solution. Therefore the Schwarz waveform Relaxation Algorithm is well defined. If it converges, 
the limit in each subdomain is denoted by Vi- It satisfies the same scheme as U at initial time, in the 
interior and on the exterior boundaries. At point it satisfies for any n > 

B-(f,g-)Vi(0,n)=B-{f,g-)V i - 1 (Ji- 1 ,n) , B+(f,g + )V i - 1 (J i - U n)=B + (f,g + )V i (0,n). (44) 



Theorem 5.1 Suppose that f is affine in the third variable d x u. Then, if the discrete algorithm con- 
verges, it converges to the discrete approximation U of problem Q). 

Proof Since p is in H 2 (£l) and q is in i? 1 (J7), there are both continuous and we have for any i Vi(0, 0) = 
-Pj(O) = Vi_i(Ji_i,0) = Pi_i(Ji_i) =p(ai), and Qj_i(Jj_i) = Qj(0) = q(a,i). We write the transmission 
conditions (|4~4"|) for n = 0. The nonlinear terms containing g^ on both sides cancel out, and we have 

A qn /XT /XT 

(D+ -D+ +—Dt)V i (0,0)-—q(a i )-—f(p(a i ),q(a i ),D+p(a i )) = 

(D+ -D- -^D+)V i - 1 (J i - 1 ,0) + ^q(a i ) + ^f(p(a i ),q(a i ),D-p(a i )), (45) 

(Dt +D+ -^Dt)V i (0,0) + ^q(a i ) + ^f(p(a i ),q(a i ),D+p(a i )) = 

\ j- /\ ftp /\ np 

(D+ +D- + —Dt)V i . l {J i - l ,Q)-—q{a i )-—f{p{a i ),q{a i )),D-p{ ai )). (46) 

Adding gSJ and @BJ> yields 1^(0,0) = L> f + V*_i(J 4 _i,0), and hence ^(0,1) = We 
define now V"(j, n) for < j < J as V^n) = V^(j, n) if j'Ax = a* + jAx, with 1 < j < J;. With the 
assumption on /, since D® = (Z?+ + D~ )/2, we can rewrite (|45l) as 

(£>+ -D" -2^D+)V(j l ,0) + 2^q(a l ) + Axf(p(a i ),q(a l ),D x p(a l ))=0, 

with = Ji + . . . + Ji, which, multiplying by — proves that V is solution of ([36]) at any point, and 
therefore U and V coincide at time and 1. A simple recursion with the explicit schemes now proves 
that, for any n, for any i, Vi(Q, n) = Vi-i(Ji-x, n), and therefore U = V. ■ 



Remark 5.2 The assumption on f in Theorem 15. JI is fulfilled when f(u,u t ,u x ) = fi(u) + f2(u)ut 
h{u)u x . 
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(49) 



6 Convergence of the Discrete Algorithm 

According to Theorem 15-11 we suppose that / is affine in d x u. We introduce the linear transmission 
operators defined by 

T ± = B ± (0,0), T ± = B ± (0,0) 
We note Ui = £//sv With these notations, the error U k = U k — Ui is solution of the linear problem 

{D+ D- t - D+ D- + l)U k = F k , in [1, Jj — 1] X [1, N] (47) 

(48) 

with the initial value U k (j, 1) = U k (j, 0) = and the transmission conditions 

T-U?(0, •) + GT' k = f-u£?(Ji-i, + G-^T\ in \0,N\, 
T+U?(J U •) + G+> k = T+U&Q, •) + G+L*- 1 , in [0,7V]. 
The remainders are given by 

F? = f(U?, U k ,D° x U k ) - f(U u Dr U k ,D° x U) + U k for n > 1, 
G-- k = —ft (0, •) + g-(U?(0, •)) - 9(Ui(0, •)), 

G^ll = ^flM-xr) + 9-{Ulx{Ji-u-))-9-{Ui-M-xr)), (50) 

Gf' k = -^fiVi* •) + 9 + (Ut(^ ')) - 9 + (Ui(Ji, •)), 

Gt+l = ^/f + i(0, •) + <? + (^ +1 (0, •)) - .9 + (£W0, •)), 

and Gf' k (0) — 0, G,f' fc (0) = 0. For n = 0, the centered derivative in x are replaced in the expression of 
F k by a forward or backward derivative. 

We define now a discrete energy as follows. We consider sequences of the form V = {^(i)}o<j<J m 
R J+1 , and we define a bilinear form on by 

J .7-1 

a A (V, W) = Ax(J2 D- (V)(j) ■ D- (W)(j) + £ V(j)W(j)). (51) 

3=1 3=1 

For a mesh function V of time and space, we define 

A* 7 " 1 



E K {V)(n) = —J2i( D t V{ 3l n)f + (D< V(J, n + l)) 2 ), 

3=1 

E P (V)(n) = a A (V(;n),V(;n-l)), 
E = ~\~ Ep. 



(52) 



The quantity -E^ is a discrete kinetic energy. It is less evident to identify Ep as discrete potential energy. 
The following lemma gives a lower bound for E under a CFL condition, and hence shows that E is then 
indeed an energy. The proof is classical ([7]) and is omitted here. 



A f2 A 4-2 \ 

E(V){n) > ( 1 - — - — E K {V)(n). (53) 



Lemma 6.1 For any n > 1, we have 

n\ > I 1 - 

As 2 4 
Hence, under the CFL condition 

At 2 At 2 1 

A^ + — <1 ' (54) 

_E is bounded from below by an energy. 



14 



The following energy estimate is obtained by a discrete integration by parts: 
Lemma 6.2 For any V solution of 

(D+D7 -D+D- +1)V = F, 1<J< J -1,1 <n<N, (55) 
we have for any n > 1, 

E{V){n) E{V)(n 1) + ^[(T + V(0, n) f + (T~ F( J, n)) 2 ] 

At ( 56 ) 
= -y l(T-V{0, n)f + (T + V(J, n)f] + 2AtAx *U n)D° V(j, n), 

i=i 

and for n = 0, 

E K (U)(0) + E(U)(0) + ^[(f+U(0, 0)) 2 + (f-U(J + 1, 0)) 2 ] 

At (57) 
= _[(T-[/(0,0)) 2 + (T + U(J,0)) 2 } +a A (P,P) + 2AxY l Qti)DtUU,0). 

i=i 

These estimates are obtained by multiplying (|55|) with V(j, n) and integrating by parts [7]. 
We now state the main result of this section. 

Theorem 6.3 Suppose that f is affine with respect to d x u. Defining the quantities 

R k (n) = G> fc (r i )(G^ fc (ri)+2T+^^^ 

-G7>\n){G7> k {n) + 2T-0?(O,n)) - Gp k (n)(Gp k (n) + 2T+U k (J h n)), 

and assuming that there exists a positive constant M such that for any iteration number K , any domain 
f2j and any discrete time n, the following estimate holds 

J 1 
2Ax J2FiU,n)D° t U k (j,n) + -R k (n) < ME(U k )(n), (59) 

j=i 

then, for At sufficiently small, the discrete Schwarz algorithm converges in the energy norm. 

Proof We apply (j56T57|) to U k . Since the initial data vanish, every term in ([FT]) vanish. Thus E(U k )(0) = 
0, and we rewrite (I56I) as 



E(U k )(n) - E(U k )(n 1) + ^-{(T+U k (0,n) + G+' fe (n)) 2 + {T-U k {J t ,n) + G~< k {n)f] 

= ^l(T~U k (0,n) + G7> k (n)f + (T+U k (J t ,n) + Gp k (n)f} (6Q) 

■ h At 
+ 2AtAxY,fKh n)D° U k (j, n) + —R k (n). 

We now insert the transmission conditions (|49|) . translate the indices in the righthand side, and add the 
contributions of all subdomains. We define a total internal energy and a total boundary energy as 

E k (n) = ]T E{U k ){n), E k (n) = ^ £[(T+f7f(0, n) + G+< fe (n)) 2 + (f~U k (Ji,n) + G7' k {n)f]. 

i—l i 

With these notations we can write 

i J, 



E k (n) - E k (n 1) + E k (n) < E k B \n) + ^AtAxJ^fi U k (j,n) + ^R k (n)). (61) 

i=l j=l 
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We now sum up JBTJ for 1 < k < K, and define Ef(n) = Y,k=i E^{n): 

K I Ji 1 

Ef{n) - Ef{n 1) + E*(n) < E° B (n) + At £ £[2Ax £ /* (j, n)£>° E^.n) + ±itf(n)]. (62) 

fc = l 8=1 J=l 

Under assumption ([59]) we deduce that 

Ef(n) - Ef(n - 1) < F (n) + MAtEf(n). (63) 

The recursive inequality is easy to solve. For At sufficiently small, MAt < 1, and since Ek(Q) — 0, we 
get 

n N 

Ef{n) < e MnAt MP) < e MT ^ F (p). (64) 
p—i p—^ 

This proves that sup n<N Ef (n) is bounded as K tends to infinity. Therefore we have 

Vn<N,Y^ E(U* - U){n) -> as k -> +oo (65) 

i=l 

which concludes the proof. ■ 

We are able to prove the assumption (|5T?|) in the case of the linear transmission conditions. 

Corollary 6.4 Suppose that f is affine with respect to d x u. For At sufficiently small, there exists a time 
T such that the discrete Schwarz waveform relaxation algorithm |^3| ) with linear transmission conditions 
(i.e. g — 0) converges to the discrete approximation U of problem ^ on (0,T). 

Proof Here the remainder (n) reduces to 

itf(n)= 2A X (~f^{Q 1 n)D Q t U^0 1 n) + f^{J u n)D° t U , [(J ll n)), (66) 
and the estimate ()59j) amounts to proving that for any k < K, n < N, 

Ji 

Ax J2 fi(j, n)D° t Uf{j, n) < M^Uf^n). (67) 

3=0 

If / is globally Lipschitz in both variables, this is merely an application of the discrete Cauchy-Schwarz 
lemma. I 



7 Numerical Results 

7.1 Remarks on Overlapping versus Nonoverlapping Schwarz Waveform Re- 
laxation Algorithms and Variants 

The original Schwarz algorithm uses overlapping domains (domains tti — (0^,6^) for 1 < i < I, with 
ai < a,i + i < bi < fei+i), with an exchange of Dirichlet data on the boundary, ft is known for elliptic 
problems to converge, but the smaller the overlap, the slower the convergence jSJ. Due to the finite speed 
of propagation, it converges in a finite number of iterations, given by N = - \cT/L] where c is the wave 
speed, and L the size of the overlap (see [7] for the linear wave equation). Therefore the convergence 
can be very slow. In the linear case, using absorbing boundary conditions instead of the Dirichlet trans- 
mission, even without overlap, improves drastically the convergence, giving in one dimension a number 
of iterations equal to 2 for some T [7]. In the nonlinear case, we have theoretical convergence results on 
the linear algorithm for sufficiently small T. However, we will see that the latter performs very well on a 
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large time interval, and that the nonlinear algorithm performs even better. 

Our experiments concern the space domain il = [0,4], simulating R with linear absorbing boundary 
conditions at each boundary. The time interval is [0,2]. Q is divided in two subdomains. The initial 
value is supported in (0, 2), with p(x) — x 3 (2 — x) 3 , the initial velocity is q(x) — 3x 2 (2 — x) 2 (x — 1). This 
is a good test since the solution is supported in the first subdomain at t = and escapes in the second 
domain before the end of the computation. 

Note that the Schwarz algorithm can be viewed as a fixed point algorithm applied to the interface 
problems 

(H+,H-) i * (B+U 2 (0, •), B-UtiJt, •)) = A{H+, H^). 

In all cases, the stopping criterion in the algorithm will be on the residual res k for A. We also compute 
the exact discrete solution in f2, and measure the discrete global error E k in L°°(0, T, L 2 (Vl)). 



7.2 The Classical Overlapping Schwarz Algorithm 

In this case, Q,\ = (0, 2 + L) and SI 2 = (2,4). The stopping criterion pertains to the residual for the 
interface problem: 

res k = (\\{u k+1 - 4)(2 + L, .)||| a + ||(^ +1 - u k ){2, .)|||,) ' 

We run the computation until the residual res k is equal to zero. The theoretical minimal number of 

Ax T 

iterations for the discrete algorithm is N t h — 1 • 

We start with the nonlinear term f(u,u t ,u x ) = u 3 . Table [T] gives the number of iterations N comp 
needed to achieve convergence (i.e the error is zero), together with Nth- O n the left we choose the overlap 
equal to 8 grid points, and vary At = 1/120 and Ax = 1/100 to fulfill the CFL condition. On the right 
we fix At = 1/120 and Ax = 1/100, and vary the overlap L. 



Ax 


1/100 


1/200 


1/400 


At 


1/120 


1/240 


1/480 


comp 


29 


54 


105 


Nth 


30 


60 


120 



L = 8Ax 



overlap L 


2 Ax 


AAx 


8A2; 


16Ax 


N 

1 v comp 


108 


55 


28 


15 


N th 


121 


61 


31 


16 



Ax = 1/100, At = 1/120 



Table 1: Number of iterations to achieve convergence for the classical Schwarz algorithm 



We show on Figure [T] the convergence history of the classical Schwarz algorithm for various values 
of the mesh size and an overlap equal to eight gridpoints.We check in each case that the error E k 
vanishes together with the residual. Furthermore we can see that the error decays very slowly for many 
iterations, and reaches zero in a few iterations, independently of the mesh size (three or four in all cases). 
The behaviour is very similar to what happens for the linear wave equation : only the finite speed of 
propagation produces convergence, which takes place when the signal has left the domain. The algorithm 
behaves similarly for other nonlinearities. 



17 




Ax 


= 0.01 


Ax 


= 0.005 


Ax 


= 0.0025 



60 80 100 

Iteration 



10" 
10"' 
10"' 
10"' 















\ 










\ 

■ 








Ax = 0.01 






Ax = 0.005 






A X = 0.0025 











60 
llcralion 



Figure 1: Variation of the residual res k (left) and the error E k (right) as a function of the iteration 
number k. 

7.3 The Non linear Nonoverlapping Schwarz Algorithms 

We will see in this section that our strategy greatly improves the performances of the Classical Schwarz 
algorithm. Note that, whereas the classical algorithm only converges in the presence of an overlap, and 
the smaller the overlap, the slowlier the convergence, our algorithms are run without overlap, 
Here the residual is given by: 

res k = (j|(G+> fe+1 -G+> fc )( ai ,.)||i 2 + ||(G- fe+1 - G^ k )(a u -)\\lX' 2 . 

We start with f — u 3 . In this case the transmission operators are the same and linear, since g^ = 0. We 
have proved in Corollarv l6.4l that there exists a final time T for which the discrete algorithm is convergent. 

In the forthcoming computations, the theoretical and numerical data are the same as before. Figured] 
plots the convergence history for various mesh sizes. The computation is stopped as soon as the residual 
reaches 0.5 10 -7 . 




Figure 2: / = u 3 . Left residual, right: error. *: Ax = 1/100, At = 1/120, o: Ax = 1/200, At = 1/240, 
0: Ax = 1/200, At = 1/240. 

We see that the algorithm converges very rapidly, independently of the mesh size. 

We consider now the case / = u 2 u x , and nonlinear transmission conditions, i.e. g < (u) = ±it 3 /6. In 
Figure [3j we plot the convergence history for various mesh sizes. The computation is stopped as soon as 
the residual reaches 0.5 10~ 7 . 
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Ax = 1/100, At = 1/120 Ax = 1/200, At = 1/240 Ax = 1/400, At = 1/480 



Figure 3: / = u 2 u x . Top: residual, bottom: error. Solid: linear, dash: nonlinear. 

Both the linear and the nonlinear transmission conditions behaves very well. The convergence takes place 
in 5 iterations with the linear transmission condition, in 4 iterations with the nonlinear one. 

We also can vary the nonlincarities in the transmission conditions, we use a real parameter 8, and 
g ± = ±Su 3 . The nonlinear strategy corresponds to 8 = 1/6, whereas the linear one is obtained for 5 = 0. 
We draw in Figure 2] the error curves after 3 iterations for the same initial values as before, the mesh 
sizes are Ax — 1/100, At = 1/120. We observe that the nonlinear strategy corresponds precisely to the 
optimal numerical value of the parameter 8, validating the high frequency approach. 




Figure 4: / = u 2 u x - Variations of the error as a function of the nonlinearity coefficient 8 

We have carried out the same computations in the case / = u 2 u t . We do not display the results here 
since they are very similar. 

8 Conclusion 

We have presented a linear and a nonlinear Schwarz waveform relaxation algorithm without overlap 
for the semilinear wave equation. On the continuous level, we proved the convergence for sufficiently 
small time intervals. We designed a discrete algorithm, in such a way that, if convergent, the algorithm 
converges to the discrete solution in the whole domain, which we prove when the nonlinearity is affine in 
d x u. In that case we proved the convergence for the linear transmission condition. Numerical experiments 
highlight the fast convergence to the discrete full domain solution in a large time domain in both linear 
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and nonlinear strategies, without overlap. Furthermore, we have shown that our nonlinear transmission 
conditions give optimal results within a large class of transmission conditions. 
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